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ABSTRACT 


A classical control system with a full-order observer is 
designed to stabilize the motions of towed vessels. The 
control method is transverse movement of the towline 
attachment point on the towed vessel. The linearized sway and 
yaw equations of motion are developed, leading to the control 
system design. The control system is tested using MATRIX x . 
Results for a barge, a mariner-class ship and a tanker are 
presented. Possible benefits of the implementation of such 
a system include improved fuel economy, a wider range of 
environmental conditions during which towing operations can 
be conducted, and improved safety. ) ■ !j('o } 
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I. INTRODUCTION 


Towing operations sometimes result in loss of property 
and/or life, owing to the inherently dangerous nature of the 
towing operation. The instability of the tow can lead to 
uncontrolled motion of the towed and towing vessels. Much 
research has been performed concerning towed vessel 
stability, and stable configurations of the towline 
attachment point, towline length, towline tension and other 
important parameters [Ref. 1,2, and 3]. 

This paper investigates active towing control, as can be 
provided by a classical control system. Athwartship motion 
of the towline attachment point on the towed vessel can 
provide some measure of control over the motions of the tow. 
A small initial offset from the path is the "disturbance". 
Since this analysis involves linearization of the equations 
of motion, and uses only proportional control, response to a 
constant or time-variant disturbance such as wind, wave or 
current action, was not studied. Further studies utilizing 
the nonlinear equations of motion and/or proportional- 
integral control, would provide assistance in determining 
the feasibility of the installation of such a system. 
Possible benefits of the implementation of such a system 
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include improved safety, increased efficiency, and a wider 
range of environmental conditions during which towing 
operations can take place. 

An initial discussion of the development of the sway and 
yaw equations of motion is pursued, and then these equations 
are applied to the towing situation. The controller design 
is explained, and since all state variables cannot be 
measured, a full-order observer is designed. Finally, the 
system is implemented analytically using Matrix x . 

Matrix x provides the tools necessary for investigating 
the response of the single-input single-output (SISO) system 
developed in this study. The ability to implement the 
control system graphically, in block diagram form, allows 
the user to more easily and understandably make alterations 
to the control system. The on-line graphics capabilities of 
Matrix x allow almost immediate feedback concerning the 
performance. 
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II. PROBLEM DEVELOPMENT 


A. EQUATIONS OF MOTION 


The equations of motion will be developed using two 
coordinate systems. The origin of the first is fixed in 
space, while the origin of the second is at the center of 
gravity of the towed vessel. In general, the vessel has six 
degrees of freedom of motion: 

* surge along the x^-axis 

* roll about the x -axis 

g 

* sway along the y -axis 

y 

* pitching about the y -axis 

y 

* heave along the z -axis 

g 

* yaw about the z -axis 

g 

This paper addresses only motions in the horizontal plane, and 
assumes a constant surge velocity component, u. Other 
assumptions include: 

* no towed vessel/towing vessel interaction 

* massless, inextensible towline 

* small motions about the initial conditions 

* no control time-lag 

Therefore, only the sway and yaw equations of motion will be 
developed. Using Figure 1, it can be shown that 


EXg = ZX^cosf + SY^sin^r and 

(1A) 

EYg = -EX f sint + ZY f cosijr, 

(IB) 
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(1C) 


2N g= = ZN f 

where X„ are forces in the x direction and Y are forces in 

g g g 

the y g direction. Expressions for and Y f must be found. 

N g are the moments about the center of gravity. 

If a similar velocity diagram is drawn, it can be shown th*t 
x^ = ucosiji - vsinijf and 
y^ = usinijf + vcosijf. 

Differentiating with respect to time yields 
x f = -u\jrsin\|f + ucosijr - lj/vcosiji - vsini|i and 
y f = ijmcosi|f usiniji - ijrvsin^ + vcosijr. 

Or, after collecting like terms, 


x f = ucosijf - vsiniji - 4r(usini|r + vcosijr) and (2A) 

y f = usini^ + vcosi^ + ij»(ucos^ - vsin\|i) . (2B) 

From Newton's laws of motion it is apparent that 
ZX f = mx^ and (3A) 

ZY f = my f . (3B) 


Substituting equations (2A) and (2B) into equations (3A) and 
(3B) yields the desired expressions for ZX f and EY^. 
Combining these with equations (1A) and (IB) aives 
SX g = m(ucosi|r - vsiniji - ijr(usini|r + vcosf)) cosijr 

+ m(usinijr + vcosiji + ^(ucosiji - vsintjr) ) siniji 

• 2*2 • 2 • 2 
= m((ucos if - i^vcos f) + (usm ^ - \Jivsin if)) 

ZX g = m[u - vf]. (4A) 

Likewise, 
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EYg = -m(ucosi|r - vsinijr - i£(usini|i + vcosiji)) siniji 

+ m(usiniji + vcos^i + ifr(ucosi|r - vsin\)r) ) cosf 

• • 2 • . 2 *2 * 2 

= m(vsin tjt + ijrusm ijr + vcos f + \|rucos |) 

EY g = m[v + ijru] . (4B) 

From equation (1C), 

2N g - I Z V (4C) 

Equations (4A), (4B), and (4C) give the forces acting 

on the center of gravity of the towed vessel in terms of the 

motion of the center of gravity. The forces implicity 

stated in X , Y , and N are not explicity stated on the 
y y y 

right hand side of the equations. Forces such as towline 

force, propellor thrust, drag, and rudder do not appear, but 

are implicit in X, Y, and N. (For simplicity, X , Y , and N 

y y y 

will be denoted as X, Y, and N throughout the remainder of 
the paper.) Many of these forces are nonlinear, but by 
expanding force terms using a Taylor series and keeping just 
the linear terms, expressions for small deviations from the 
initial states can be written. It is assumed that the force 
terms X, Y, and N are functions of the velocities and 
accelerations u, v, u, v, if, and V • Then, 

Y = (<3Y/0u) o au +(0Y/dv) Q AV +(dY/<3u) Q AU +(dY/8v) Q AV 

+ (8Y/df) o Af +(dY/dV* ) o a? +Y q 
F rom symmetry of the towed vessel about its longitudinal 
axis it apparent that, 

OY/au) o = (aY/au) o = o. 
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So, 


Y = (3y/3v) q av + (3y/3v) o av + (3y/3^) q ai^ 

+ (3y/3*) 0 aV + V (5A) 

To obtain an equivalent expression for Y, substitute the 

perturbed values into the right hand side of equation (4B). 

Y = m[(v o + av) + (+ o + A$)(U Q + AU)] 

= [V q + AV + ^ Q U 0 + ^ q AU + A ^ U Q + A \|TA Ul ] m. 

Dropping the higher order terms yields, 

Y = m[V o + 4» 0 U q + av + i£ q au + Atjiu o ]. (5B) 

After noting that the first two terms in equation (5B) are 
Y q , equations (5A) and (5B) can be equated to show that 

(3Y/3v) q av + (3y/3v) q av + (3y/3^) o a$ + (3y/3|) o aV* = 

m[AV + ib au + a iu ] . 

Adopting the notation (3Y/3v) q = Y v , (3y/3v) q = Y^, 

(3y/3$) q = Y* , (3y/ 3\j») q = Y- , av = v, av = v, Aijr = ljr, and 

A*i|i = V/ and assuming v q = u q = v q = ij» o = $ = 0, this equation 
can be rewritten as 

Y v V + Y^ V + Yj * + Yj * = m[ v + iu o ]. (6A) 

It must be remembered that v, v, ijr, and denote small 
deviations from the initial state. Likewise, an equation 
for the yaw forces can be developed yielding, 

N v v + v + f + N” V = I z V. (6B) 

Equations (6A) and (6B) can be rewritten as 

Y v v + (Y^ - m)v + (Y^ - mu Q )* + Y-f - F a = 0 (7A) 

N v v + + (N—I z )V - M a = 0, (?B) 

where F are applied forces and M are applied moments. 
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DEVELOPMENT OF THE 2-DIMENSIONAL TOWING STATE EQUATIONS 


For the situation depicted in Figure 2, equations (7A) 
and (7B) can be written as 

Y v v + (Y^ - m)v + (Yj - mu Q )* + Y"V 

- Tsin(y + tjr) = 0 (8A) 

N v v + N^v + + (N~ - I z ) v - Tx p sin(y + i|r) 

- Ty p cos(y + i|r) =0 (8B) 

Motion of the towed vessel will be controlled by varying y . 
The terms involving the towline tension, T, are nonlinear 
and must be linearized as follows before proceeding further. 
For small angles y and i|r, 


siny = y = OP/L = (y + x i|r + y )/L and 

P P 

sin ( y + i|r) = y + i|r = (y + x p i(r + y p )/L + i|r 
= *(1 + Xp/L) + y/L + y p /L. 

Also, 

y = v + u o sini|r a v + u q i)i and 

v = y - u o * so, v = y - d(u o f)/dt = y - u q $. 
Substituting these expressions into equations (8A) and (8B) 
yields, 

(Y^ - m)y + Y v y - Ty/L + Y-* - (Y^u - Y^)$r 
-[Y y u + T (1 + Xp/L )]♦ = Ty^L and (9A) 

N^y + N y y - Tx p y/L + (N- - 1.,)+ - (N^u - 
- [N v U + Tx p (1 + X p /L)]f = TX p y p /L + Ty p . (9B) 

Then, by assigning the state variables 
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Figure 2. Problem Geometry 
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X 1 = y, X 2 = f, X 3 = y , and X = 
and substituting into equations (9A) and (9B) 

(** - m)X 3 + Y»X 4 = -Y v X 3 + TX X /L + (Y^u - Y^)X 4 + 

[Y v u + T(1 + x p /L)]X 2 + Ty p /L (10A) 

V 3 + (N * " J z )X 4 = -N v X 3 + Tx p X l /L + (N v U " % )X 4 + 

[N v u + Tx p (1 + x p /L)]X 2 + Tx p y p /L + Ty p (10B) 

From the original equations defining the state variables, 

it can be seen that 

X. = X 3 and (HA) 

X 2 = X 4 . (11B) 

These comprise the first two of four state equations. The 
second two state equations must be found from equations 

(IOA) and (10B). Multiplying equation (10A) by (N” - l z ) 
and equation (10B) by (-Y”) and adding yields an equation 
for X 3 . 

[(%* - V (Y v " m) - ¥'^3 = 

T/L[(N” - I z ) - X p Y-]X 1 + 

{[Y v u + T(1 + Xp/L)](N^-I z ) - Y»[N y u + Tx p (l + x p /L)]}X 2 
+ [N v Y» - Y y (N^ - I z )]X 3 
+ C(V - Y-)(Nf - I z ) - Yy(N^u - N^)]X 4 
+ T((N- - I z )/L - Y»(l + X p /L)]y p . (11C) 

Similarly, multiplying equation (10A) by (-N^) and equation 

(IOB) by (Y^ - m) and adding yields an expression for X 4> 
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C( N * - V - ■) - N * Y $]X 4 = 

T/L[X p (Y^ - m) - N <r ]X 1 + 

{(Y^ - m)[N v U + Tx p (1 + X p /L)]-N^[Y v U + T(1 + X p /L)]}X 2 

+ [N v Y v " N v< Y * " m M x 3 
+ (-N^mu - N # Y^ + N^m + N^Y-)X 4 

+ T[(Y^ - m)(l + x p /L) - N^/L]y p (11D) 

Equations (11A), (11B), (11C) and (11D) can be more 
conveniently expressed in the matrix form 



Where, 

D = (Ny - X z ) (Y^. - m) - 
a 31 = d'VlKNj - I z ) - Xp Y ¥ ] 

a 32~ D' 1 ([Y V U+T(l+X p /L)] ( N ¥ - I z , -Y ¥ [H V U + Tx p (l +x p /L)]} 

a 3 tr D ‘ 1 IV» - Y v< N r I z )1 

a 34 = D _ 1 [(Y^U - Y.)(N r I 2 )-Y ¥ (N <> u - H + ) ] 

a 4i = D _ 1 T/L[x p (Y 4 - n)-M^] 

a 42 = D" 1 {(Y i> -m)[N v u + Tx p (l+x p /L)]-N^[Y v u + T(l+X p /L) ] ) 

a 43 = D_1 tV V ' V Y * - m) ) 

a 44 = D _ 1 [N^ - NjY^ + (N-- N^u)m] 
b 3 = D' 1 T[(Nj-I 2 )/L - Y^fl+Xp/L)] 
b 4 = D' 1 T[(Y^ - m)(l+x p /L) - N^/L] 
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Figure 3 is a block diagram representing the physical 
system. For implementation with Matrix x , each element of 
the A and B matrices represents a gain block. 
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Figure 3. riant Block Diagram 
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C. CONTROLLER DESIGN 


If the control, y , is proportional to the state of the 
physical system, {X}, then y p can be expressed as 
y p = ~{K}{X} 

= -K-X- - K_X_ - K_X_ - K.X. 

11 2 2 31 44 

where {K} are the proportional gains. Then, the system of 
equations represented by equation (12) can be rewritten as 
{X} = [A]{X) - {B}{K}{X} (13) 

= ([A] - (B}{K}){X) (13A) 

Taking the Laplace transform of equation (13) yields, 
s{X(s)} = [A](X(s)) - (B){K){X(s)}, 
or, 

s{X(s)} - [A]{X(s)) + {B){K)(X(s)) = {0} 

So, 

(S[I] - [A] + {B)(K})(X(s)) = (0) and 
([A] - {B}{K} - s[I])X(s) = (0), 
where [I] is the identity matrix. The closed-loop 
eigenvalues of this system are given by 

|[A] - (B)(K) - s[I]| = 0. (14) 

Since the pair (A,B) is, in general, controllable, 
values of K can be chosen which will yield any desired 
eigenvalues. Applying equation (14) to the towing equations 
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of motion yields, 
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If the desired eigenvalues are s^, s 2 , s^, and s 4 , then 
the desired characteristic equation is 


(s-s^ (s-s 2 ) (s-s 3 ) (s-s 4 )=0. (15) 

For a stable system, one in which the poles of the transfer 
function are in the left half of the complex plane, the 
eigenvalues must have negative real parts. Values of K can 
be found by equating the s-term coefficients of equation 
(14) with those of equation (15). The control law can then 
be written as 



k l k 2 k 3 k 4> 



(16) 


D. OBSERVER DESIGN 

The control, y , is a linear function of the state 

P 

variables X 1 , X 2 , X 3 , and X 4 « These state variables 
represent the cross track error from the straightline path, 
the yaw angle, the sway velocity and the yaw angular 


15 




velocity. Though required as input to the controller, one 
or more variables may not be measurable. In particular, it 
may be difficult to measure the y-distance from the path and 
sway velocity. Both yaw angle and yaw angular velocity 
could be measured, but for the purposes of this research, it 
is assumed that only the measure of the yaw angle is 
available. In this case, estimates of the remaining three 
state variables must be made and provided as input to the 
controller. Since three of the four variables are unmeas¬ 
urable, a full-order observer will be designed. The 
measurable output, y, of the physical system (plant) can be 
described by 

(y> = {C}{X}. (17) 

Upon expanding this system of equations 

* - < C l C 2 C 3 C 4> < X 1 X 2 X 3 X 4 |T ‘ 

=(010 0 }(X x x 2 x 3 x 4 ) t 

If the approximations of X^, X 2 , X 3 , and X 4 are denoted by 

AAA A 

X 2 , X 3 , and X 4 , then the control law becomes 
y p = ~{K){X}, (18) 

and the plant model can be described by 

A A 

(X) = [A]{X} + {B}y p . 

The plant model can be corrected by feeding back the 
difference between the actual plant output and the plant 
model output. The result is the observer equation, 

(X) = [A]{X} + {B)y p + {L)({y} - {C}{X}), (19) 
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where {L} is the observer gain vector. The error between 
the plant and the plant model is 

{X}-{X}=[A]{X}+(B}y p -[A]{X}-{B}y p -{L}({y}-{C){X}) 

=[A]({X}-{X})-{L){C}({X)-{X}) 

=([A]-{L}{C})({X}-{X>). 

It is then clear that the error dynamics are governed by the 
eigenvalues of 
[A]-<L}<C}. 

It is also desirable that the error die out quickly, with 
respect to the dynamics of the plant. Since 

|[A]-{L){C}-s[I]|=0, (20) 

it is possible to pick values of L such that the resulting 
eigenvalues fall to the left of the plant eigenvalues in the 
s-plane. If it is noted that the eigenvalues of 
[A] ~ { L} {C) 

are the same as the eigenvalues of 

([A]-{L){C)) T =[A] T -{C} T {L} T , 
then equation (20) can be rewritten as 

|[A] T -{C} T (L) T -s[I]|=0, (21) 

and values of L can be determined using the same method as 

that used to determine the values of the controller gain, 

{K). The block diagram of the observer is shown in 

Figure 4, and a block diagram of the complete system 

(plant, observer, controller) is shown in Figure 5. 
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Figure 4 
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Observer Block Diagram 
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Figure 5. Plant/Compensator Block Diagram 
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E. MATRIX^ IMPLEMENTATION 

The implementation of a single-input single-output 
system using Matrix x is straight forward. The coefficients 
in the [A] matrix, and the {B}, {C}, {K} and {L} vectors are 
constant with respect to time. The [A] matrix is determined 
by the hydrodynamic derivatives associated with a particular 
vessel, the towline tension, the towline length, the 
distance from the center of gravity forward to the towline 
attachment point, and the velocity component, u. The (C) 
vector is simply determined by which outputs are measurable. 
The {K> and {L} vectors are determined using the "poleplace" 
command available in Matrix x . The user can input the closed 
loop observer and controller poles interactively. Quite 
simply, Matrix x takes the drudgery out of the linear algebra 
associated with developing a control system. 

Of more interest perhaps, are the System-Build and 
Simulate functions of Matrix x> The block diagrams im¬ 
plemented using Matrix x have been included in the Appendix. 
The initial condition of the integrator between X 1 and 
was set to 0.2, simulating an off-path condition. The 
effect of the control system was to bring the vessel back on 
path. Plots of y, v, y^and are immediately available as 
output. 

Program listings for the control system with observer 
and the control system without the observer are included in 
the Appendix. 
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III. RESULTS 


The results are presented graphically, and with the 
exception of yaw angle, are in nondimensional form. Yaw 
angle is in degrees. The relationships describing the 
nondimensional quantities are included in the Appendix. 

Three vessels were studied: 

1. 191.6 foot barge 

2. 1066.3 foot tanker 

3. 528 foot mariner 

A summary of the hull particulars and hydrodynamic 
derivatives for each of these vessels is included in the 
Appendix. The mariner and the barge are stable at the four 
knot forward velocity used in this study. The tanker is 
unstable. A nondimensional towline length of 2.5 was used 
throughout. The nondimensional towline tension was taken 
from resistance curves for each of the vessels at four 
knots. The nondimensional longitudinal distance from the 
center of gravity to the towline attachment point was 
assumed to be 0.5, and is also constant throughout. The 
observer poles were selected to be five times more negative 
than the controller poles. Obviously, the towline 
attachment point can move only a finite distance to port and 
starboard. This distance is ultimately limited by the beam 
of the vessel in question. Various values were used, but 
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the maximum value used was 0.1. Reducing this control 
saturation point sufficiently will result in the inability 
of the controller to stabilize the response. Four plots 
were generated for each simulation. The y-distance from the 
path, the towline attachment point lateral offset, the sway 
velocity and the yaw angle are all plotted versus 
nondimensional time. Where the observer is implemented, the 
estimates of y, v, and ljr are incorporated in the graphs of 
y, v, and \|r. 

A. MARINER 

1. Figures 6 and 7: With Observer/Without Control 

The block diagram used in this simulation is shown 
in Figure 6. Towline attachment point movement is limited 
to zero in the controller. These graphs illustrate the 
straight line stability of the mariner, with no control. 
Notice that y, v, and ljr all oscillate about zero before a 
straightline path is achieved. The straightline path 
condition is reached in a nondimensional time of about 60. 

2. Figures 8 and 9: Without Observer/With Control 

Effectiveness of the controller is directly 

related to the accuracy of the state estimates. The block 
diagram used in this simulation is shown in Figure 8. The 
plant outputs are all assumed measurable, and are used as 
inputs to the controller. Notice that there is no overshoot 
in y, v, and f, and a straightline path is achieved in a 
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Figure 9. Mariner Without Observer/With Control 












nondimensional time of about 25. The maximum lateral offset 
of the towline attachment commanded by the controller is 
about 0.014. 

3. Figure 10: With Observer/With Control 

The block diagram used for this simulation is 
Figure 6. However, lateral towline attachment point 
movement is now limited to -0.1 to 0.1. A straightline path 
condition is reached in a nondimensional time of 20 to 25, 
but almost two times as much control effort is required as 
was without the observer. Also, note that there is some 
small overshoot in y, v, and f. Clearly, initial errors in 
the state estimates have degraded the performance of the 
system. 

B. BARGE 

1. Figure 11: With Observer/Without Control 

The block diagram used for this simulation is 
Figure 6. Again, lateral movement of the towline attachment 
point is limited to zero. Although the barge has 
straightline stability, note the errors in the estimates of 
y and v which are apparent in the graphs. A straight line 
path is achieved in a nondimensional time of about 20, but 
considerable oscillation in v and f takes place. Also of 
note is the fact that the yaw angle reaches a maximum value 
of about six degrees. This is significant. There is small 
overshoot in v and f. 
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Figure 12: Without Observer/With Control 

The block diagram used in this simulation is 
Figure 8. All plant outputs are again assumed to be 
measurable and applied as inputs to the controller. A 
straightline path is achieved in a nondimensional time of 
about 25. The maximum yaw angle has been reduced to about 1 
degree, at the expense of considerably more overshoot in y, 
v, and if. Note that the maximum lateral offset of the 
towline attachment point is about 0.029. 

3. Figure 13: With Observer/With Control 

With the observer, about half as much control 
effort is required, but the maximum yaw angle reaches 
approximately 4.5 degrees. Overshoot in y and v is 
comparable to the simulation without the observer, and a 
staightline path is attained in a nondimensional time of 
about 25. 

4. Figure 14: With Observer/With Control 

This figure illustrates the effect of reducing the 
maximum allowable control effort to less than that which was 
used in Figure 13. Since the barge is stable, the control 
does not saturate at the 0.1 limit. For this simulation, 
the maximum lateral offset of the towline attachment point 
was reduced to 0.0045. A straightline path is achieved in 
about the same time as before, but the maximum yaw angle 
increases to about 5.5 degrees and the control is saturated 
twice. Higher sway velocity also results. 
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C. TANKER 

Since the tanker is unstable, it makes for the best 
test of the control system. A test of the robustness of the 
control system is also performed, both with and without the 
observer. 

1. Figure 15: With Observer/Without Control 

This figure illustrates the unstable response of 
the tanker in the absence of control. 

2. Figure 16: Without Observer/With Control 

The control system achieves a straightline path in 
a nondimensional time of about 25, with no overshoot in y, 
v, or if. The maximum yaw angle attained is about 1.4 5 
degrees and the towline attachment point only moves 0.06 
ship lengths laterally. However, 0.06 ship lengths in 
dimensional form is approximately 64 feet for the tanker. 
This much lateral movement would reguire a beam in excess of 
126 feet. 

3. Figure 17: With Observer/With Control 

For this simulation, a maximum towline attachment 
point lateral offset of 0.1 was used. Again, a straightline 
path was attained in a nondimensional time of about 25, but 
the yaw angle exceeded three degrees and the control was 
saturated twice. There is also overshoot in y, v, and 
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Figure 16. Tanker Without Observer/With Control 
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Figure 17. Tanker With Observer/With Control 
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4. Figure 18: With Observer/With Control 

The control saturation limit was set to 0.05 for 
this simulation. A beam of greater that 107 feet would be 
required for the traversing mechanism. It takes 
considerably longer to attain a straightline path and the 
control is saturated three times. The overshoot in y, v, 
and i|r is greater. However, the maximum yaw angle attained 
is less than was observed with a saturation limit of 0.1. 

5. Figure 19: With Observer/With Control 

This simulation illustrates the effect of reducing 
the control saturation limit. The limit was set to 0.035, 
corresponding to a required beam of about 75 feet. Note 
that the response is unstable. In defense of the control 
system, it might be noted that the initial offset from the 
path of 0.2 ship lengths (213 feet) is significant. 

6. Figure 20: Test For Robustness 

This simulation was performed with the observer 
implemented and with control. It was assumed that the value 
of dY/dv was in error. This hydrodynamic derivative 
represents damping. Reducing the damping in the y-direction 
should destabilize the response to some extent. The 
original controller and observer gains were retained, but 
for purposes of modelling the plant, the reduced value was 
used. The magnitude of dY/dv was reduced by ten per cent. 
Clearly, the system is less stable than before, and takes a 
nondimensional time of about 50 to achieve a straightline 
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Figure 18. Tanker With Observer/With Control (0.05) 
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Figure 19. Tanker With Observer/With Control (0.035) 
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path. The overshoot in y, v, and f is significantly 
greater, and a maximum yaw angle of almost five degrees is 
attained. 

7. Figure 21: Test For Robustness 

As mentioned previously, errors in the state 
estimates adversely impact the ability of the controller to 
stabilize the system. For this simulation the block diagram 
in Figure 8 was used. The plant outputs are assumed to be 
measurable. With a 65 per cent reduction in dY/dv, the 
controller was still able to stabilize the response easily. 
Obviously, state estimate errors are degrading the 
performance of the system. 

8. Figures 22 and 23: Without Observer/With Modified 

Controller 

For the simulations that produced the results in 
these figures, the controller was modified. The gains 
acting on y and v were set to zero in the controller. This 
represents the inability to measure y and v. Then, it was 
assumed that both f and ijr could be measured, and applied as 
inputs to the controller. The state variables y and v were 
not estimated, so y p was dependent solely on ijr and f. 

Compare Figures 16 and 22. The difference between these two 
simulations is that in the first, y p is dependent on y, v, 
t, and f. In other words, all four state variables are 
assumed measurable and applied as inputs to the controller. 
In the latter, y p is dependent on only f and f. It is 
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apparent that the controller performs in an almost identical 
fashion without v and y as inputs to the controller. 

Compare Figures 18 and 23. For both, the control is 
saturated at 0.035. The simulation represented by Figure 18 
used the full order observer, assuming only if was 
measurable, and providing estimates for if, y, and v. The 
performance of the modified controller, without estimates or 
measurements of v and y, is obviously much faster, with 
almost no overshoot in if, y, and v. It seems clear then, 
that not only does estimation of y and v fail to improve the 
performance of the compensator, it actually degrades its 
performance as compared to a system assuming measurement of 








IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This study investigated the possibility of athwartships 
movement of the towline attachment point as a means of 
improving the stability of towed vessels. The assumptions 
with respect to model linearity restrict its applicability 
to trajectories close to the path of the towing ship. 
Additionally, the controller design in combination with the 
full order observer did not prove to be exceptionally 
robust. In the aosence of the observer, and assuming that 
measurements of i|i and ijr were possible, the compensator 
performed admirably and was able to tolerate a 65 per cent 
reduction in the value of the damping hydrodynamic 
derivative, Y v . Determining whether or not lateral movement 
of the towline attachment point would provide a sufficient 
measure of control in and of itself, will require further 
study. It has been established in the literature that there 
are an infinite number of combinations of x p and towline 
length that will provide for stability of the tow. 
Unfortunately, not all combinations of x p and towline length 
are feasible. For the tanker, an x p of greater than one 
would be required to ensure stability of the system with a 
nondimensional towline length of 2.5. This study has shown 
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that for the tanker, the distance y p required to provide an 
acceptable measure of improvement in stability may be too 
large to be accommodated by the beam of the vessel. This 
distance may actually be larger than predicted, because no 
delay between command and movement of the attachment point 
was modeled. Some combination of increased x and movement 

p - 

of the towline attachment point would probably be a more 
effective means of control than either individually. 
Practical considerations may preclude the actual 
implementation of such a system. 

The performance and robustness of the system was 
seriously degraded by errors in the state estimates. This 
was most evident in the barge simulations, and in the tests 
for robustness performed with the tanker model. The data 
indicate that the best compensator performance is obtained 
in the absence of state estimates for v and y. Measurement 
of and ^ is possible, and in fact, the presence of 
measurements of v and y provided little improvement over 
simply measuring iji and ij». Should further research prove 
this to be so, the cost and complexity associated with the 
entire observer can be avoided. 

B. RECOMMENDATIONS 

Modelling of the time delay between control command and 
response would result in a further increase in the amount of 
control effort required. As mentioned previously, some 
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combination of an increase in x p and lateral movement of the 
towline attachment point may be workable. At any rate, an 
analysis based on the nonlinear equations of motion and 
modelling of the control time delay, would help to determine 
the feasibility of the system. The most apparent 
disadvantage of a y p -only compensator, would seem to be the 
sheer size of the traversing mechanism required to provided 
the necessary control effort. Any reduction in the amount 
of control effort required, whether through increased 
towline tension or an increase in x p , would improve the 
performance of the control system. 

Additionally, to provide true pathkeeping stability in 
the presence of constant or time-variant disturbances, 
integral control needs to be implemented. Krikelis and 
Kavouras showed that the path of a towed vessel could be 
maintained within a narrow band using a nonlinear rudder 
control, in the presence of oblique currents of one knot, 
and wave action.[Ref. 4] Given the effectiveness of a 
rudder at low speeds, it should be possible to provide an 
equal level of control through movement of the towline 
attachment point. 

This study has shown that it is possible to stabilize 
the motions of a tow with lateral movement of the towline 
attachment point on the tow. Further research is necessary 
to better model the towed vessel dynamics and the 
command/actuator time delay. Nonlinear analysis of the 
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towed vessel dynamics will shed more light on the 
feasibility of the system. Integral control may be 
implemented to provide true pathkeeping ability, although 
simply stabilizing the response will improve safety. 
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//****(,*.**** ft cmacinerwol . x a*************************,* *// 


// // 

// Matrixx variables for use in designing // 

// mariner towing classical control system. // 

// // 

//The following are the non-dimensional hydrodynamic // 

//coefficients, unless otherwise indicated. // 

// // 


Xudot“-0.000444; 

M = 0.00888 ; 

Yvdot--0.00912; 

Y v « - 0.0 1 4 34 ; - 
Yrdot-0.0; 

Yr~0.00456 ; 

Yde1=0.0278; 

Nvdot»0.0; 

UV--0.0046; 

Nrdotlz=-0.00115; //Nrdot-Iz// 
Nr--0.00296 ; - 
Nde1--0.0139; 


// // 

//End list of hydrodynamic derivatives // 

// // 

// // 

//The following statements define the A and B matrix elements.// 

// // 

//***“** Get critical parameters from user ***»*****»»»***// 

// // 

inquire LL //nondimensiona1 towline length // 

inquire Xp //nondimensional distance from CG fwd to towline...// 

//attachment// 

inquire T //nondime nsiona1 towline tension // 

inquire YpSat //nondimensional max control// 


// // 

U=1.0; //nondimensional ship velocity// 

L-1.0; //nondimensional ship length// 

// // 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

// // 

//****»*»*.*«**+**** A matrix elements *.****++*»*******»*// 

// // 

a 13-1.0; 


a 2 4 = 1 . 0 ; 


a 31-( 1/D) MT/LL) *(NrdotIz-Xp*Yrdot) ; 

a 3 2 = (1/D)* ( (Yv *U + T*(1+Xp/LL))*NrdotIz-Yrdot*(Nv*U+T*Xp*. . . 
(1+Xp/LL))); 

a 3 3 = (1/D)*(Nv*Yrdot-Yv*NrdotIz) ; 

a 3 4 -(1/D)*( (Yvdot *U-Yr)‘NrdotIz-Yrdot*(Nvdot *U-Nr) ); 
a 4l-(1/D)*(T/LL)*(Xp*(Yvdot-M)-Nvdot); 
a 4 2=( 1/D)*((Yvdot-M)*(Nv*U+T*Xp*(1 + Xp/LL))-Nvdot* . . . 

!YV *U + T*( 1+Xp/LL) ) ) ; 

a43=(l/D)*(Nvdot*Yv-Nv*(Yvdot-M)); 

a 4 4 = (1/D)*(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U)*M); 


// // 

//+*+*+******** b matrix elements **<,« + **»** + »***»*«***.»» - 

// // 

b3-(1/D)*T*(NrdotIz/LL-Yrdot*(1+Xp/LL)); 
h4-(l/D)*T*((Yvdot-M ) *( 1+Xp/LL)-Nvdot/LL) ; 

// // 

//**»+****** Matrix operations: design a controller »*+**++// 


Figure 35. Matrix^ Program Listing With Observer 


63 





// // 

A-(0,0,a 13,0;0,0,0,a24;a31,a32,a33,a34; . . . 
a41,a42,a43,a44); 

// // 

//**** Check for stability of uncontrolled system •*******»// 

//**** (Open loop eigenvalues) *************************** // 

// // 

AEVALUES-eig(A); 

// // 

//**************** Define other matrices ******************// 

// // 

B - [ 0.0 ; 0.0 ; b 3 ; b 4 ] ; 

C»|0,1,0,0); //Note: assuming only yaw angle available// 

DMa trix-10 ) ; 

SysMat-[A,B;C,DMatrix ) ; //System matrix: Xdot-AX+BU, Y-CX+DU// 

// // 

//* *************** Check for controllability **************// 

// // 

NS-4; //Number of states// 

(SC,NSC,Tl)-cntriable!SysMat,NS); 

// // 

//SC-system matrix of controllable subsystem// 

//NSC-number of states controllable// 

//Tl-ttansformation between original and new space// 

// // 

//*********» Get desired closed-loop poles from user *»*»»*// 

// // 

inquire Poles 

// // 

//************ Pole place to find the K row vector ******•»// 

// // 

K = poleplace(A,B,Poles) ; 

// // 

// // 

K11 = K( 1,1); 

::12-K(1 ,2); 

K13-K(1,3); 

K 1 4 - K ( 1,4 ) ; 

// // 

//************** Check for observability *****************// 
// // 
(SOBS,NSOBS,T2)-cbservable(SysMat,NS) ; 

// // 
//SORS=observable subsystem matrix// 

//NSOBS=-numbe r of states observable// 

//T2«transformation between original and new space// 

// // 
//****»•**»** Design full order observer *****************// 
// // 
//** observer equation: Xhatdot-A*Xhat+B*u+L*(y-C*Xhat) **// 

// // 

//*********** ( 3 e (- observer poles from user *»*****»***»»*»// 

// // 
inquire obspoles 
LT-poleplace(A',C'.OBSPOLES); 

// // 

//LT-L transpose// 

//OBSPOLES-observer poles// 

// // 

//»»****** Define L matrix elements **********************// 
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// // 

111—LT(1,1); 

121* LT(1,2); 

131-LTll, 3); 

141*LT(1,4); 

// // 

//*** Define gains for use in simulation *****************// 

// // 

G2 = a32-131 ; 

G6«*a42-141; 

// // 

// *************** SIMULATION *************************** // 

// // 

sim('analyze/csystemwo/'); 

// // 

inquire tsim //simulation duration in nondimensiona1 time// 
inquire tstep//simulation time step, nondimensiona1// 

// // 

TS*(0:tstep:tsim)'; //the "time" vector for simulation// 
y-sim(TS); 

// // 

// *************** PLOTTING ****************************** // 

// // 


plot(TS,V(t,t 3 7]),'upper left xlabel/T/ ylabel/Y Yhat/TITLE/... 
Poles=[-0.5,-0.51,-0.52,-0.53)/') 

plot(TS,Y(:,5),'upper right xlabel/T/ ylabe1/Yp/TITLE/... 
ObsPoles»(-2.5,-2.55,-2.60,-2,65)/') 

plot(TS,Y(:,(4 91),'lower left xlabel/T/ ylabel/v vhat... 
/TITLE/T=0.001 Xp-0.5 LL-2,5/') 

plot(TS,Y( : , l6 8)),'lower tight xlabel/T/ ylabel/Psi Psihat... 
/TITLE/Mariner/') 
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//.*»*+******* cmarinerl.x ************************»******// 


// // 

// Matcixx variables for use in designing // 

// mariner towing classical control system. // 

// // 

//The following are the non-dimensional hydrodynamic // 

//coefficients, unless otherwise indicated. // 

// // 


Xudot--0.000444 ; 

M-0.00888; 

Yvdot>-0.00912; 

Yv--0.01434; 

Yrdot-0.0; 

Yt-0.00456; 

Ydel-0.0278; 

Nvdo t = 0.0; 

Nv--0.0046; 

Nrdotlz=-0.00115; //Nrdot-Iz// 
Nt=-0.00296; 

Ndel--0.0139 ; 


// // 

//End list of hydrodynamic derivatives // 

// // 

// // 


//The following statements define the A and B matrix elements.// 


// // 

//*.***** Get critical parameters from user +**+*+******.*»// 

// // 

inquire LL //nondimensional towline length // 

inquire Xp //nondimensional distance from CG fwd to towline...// 

//a t tachment// 

inquire T //nondimensional towline tension // 

inquire YpSat //nondimensional max control// 

// // 

U-1.0; //nondimensional ship velocity// 

L-1.0; //nondimensional ship length// 

// // 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

// // 

//*.»*»*******+*»*** a matrix elements a******..****...***// 

// // 

al3 = 1 . 0 ; 

324-1.0; 

a 31 = ( 1/D)*(T/LL)+(NrdotIz-Xp*Yrdot ) ; 

a 3 2 = < 1/D)*((Yv*U+T*(1+Xp/LL) )*NrdotIz-Yrdot*(Nv*U + T*Xp*. . . 
(1+Xp/LL))); 

a 3 3 = (1/D)*(Nv*Yrdot-Yv*NrdotIz); 

a 3 4 — (1/D)*((Yvdot*U-Yr)*NrdotIz-Yrdot*(Nvdot‘U-Nr) ) ; 
a 41 = (1/D)*(T/LL)*(Xp*(Yvdot-M)-Nvdot); 
a 42=(1/D)*((Yvdot-M)*(Nv*U + T*XpM1 + Xp/LL) )-Nvdot*... 

(Y V * U + T*( 1 + Xp/LL) ) ) ; 

a43-(l/D)*(Nvdot + Yv-Nv*(Yvdot-M) ) ; 

a44-(1/D)*(Nvdot * Yr-Nr *Yvdot+(Nr-Nvdot *U)*M); 


// // 

//***»***+***+* B matrix elements **.*****».*****..***.**.// 

// // 

b 3 —(1/D) *T* (NrdotIz/LL-Yrdot *(1 + Xp/LL )); 
ta 4-( 1/D) *T+ ( (Yvdot-M)*( 1+Xp/LL)-Nvdot/LL); 

// // 

//***»*****+ Matrix operations: design a controller **+++**// 


Figure 36. Matrix^ Program Listing Without Observer 
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// // 

A=[0,0,al3,0;0,0,0,a24;a31,a32,a33,a34;... 

a41,a42,a43,a44); 

// // 

//**** Check for stability of uncontrolled system ****»*»**// 

//**** (Open loop eigenvalues) ********************.******// 

// // 

AEVALUES-eig(A); 

// // 

//*»»****»*******» Define other matrices ***************** 4^/ 

// // 

B= ( 0.0 ; 0.0 ; b 3 ; b4 ) ; 

C-(0,1,0,0); //Note: assuming only yaw angle available// 

DMa trix-[0i; 

SysMat-(A,B;C,DMatrix1; //System matrix: Xdot-AX+BU, Y-CX+DU// 

// // 

//**************** check for controllability **************// 

// // 

NS-4; //Number of states// 

(SC,NSC,TlJ-cntrlable(SysMat,NS); 

// // 

//SC=system matrix of controllable subsystem// 

//NSC-number of states controllable// 

//Tl=transformation between original and new space// 

// // 

//*»*****•** (; e t desired closed-loop poles from user *»****// 

// // 

inquire Poles 

// // 

//************ pole place to find the K row vector *»*»****// 

// // 

K-poleplace(A,B,Poles); 


// 

// 

Kll-K(l,l); 

K12 = K(1,2); 

K1 3 - K ( 1, 3 ) ; 

K14=K(1,4); 

// 

// *************** SIMULATION *************************** 

// 

sim('analyze/csystem/’); 

// 

inquire tsim //simulation duration in nondimensional time// 
inquire tstep//simu 1 ation time step, nondimensional// 

// 

TS-(0:tstep:tsim]'; //the "time" vector for simulation// 
y=sim(TS); 

// 

// *************** PLOTTING ****************************** 

// 

plot(TS,Y(:, 3) ,'upper left xlabel/T/ ylabe1/Y/TITLE/... 
Poles-(-0.5,-0.51,-0.52,-0.53)/’) 

plot(TS,Y(:,5),'upper right xlabel/T/ ylabel/Yp/TITLE/... 

No observer/') 

plot(TS,Y(4),'lower left xlabel/T/ ylabel/v... 

/TITLE/T-0.001 Xp-0.5 LL-2.5/') 

plot(TS,Y(:,1),'lower right xlabel/T/ ylabel/Psi... 
/TITLE/Mariner/') 


// 

// 


// 

// 

// 

// 


// 


// 

// 

// 
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TABLE 1 


TOWED VESSEL DATA 




Vessel 


Prooertv 

Barae 

Tanker 

Mariner 

LBP, ft 

191.56 

1066.3 

528 

X' 

-0.00136 

-0.0009 

-0.0004444 

u 

m' 

0.170 

0.0181 

0.00888 

Y! 

-0.01383 

-0.0171 

-0.00912 

V 

Y' 

-0.0153 

-0.0261 

-0.01434 

V 

Y i 

0 

0 

0 


0.00238 

0.00365 

0.00456 

N! 

0 

0 

0 

V 

n : 

-0.007285 

-0.0105 

-0.0046 

i' “Nv 
z ijl 

0.00188 

0.00222 

0.00115 

1 

N? 

-0.00128 

-0.0048 

-0.00296 

* 
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TABLE 2 


NONDIMENSIONAL TERMS 


m' = 

m/m r 

y ;= vw 

u' = 

U/ u o 

Y *= V m r 

v 1 = 

V/U o 

Y r= Y *= V< m r u o> 

u ' = 

U/(U^-L) 

Y^= Y£= Y. f ./(m r L) 

V ' — 

v/(u‘-L) 

N v = V^W 

t' = 

t/(L/u o ) (time) 

H = V< m r L ) 

*z = 

v<v 2 > 

N r = N i = V (m r U o L) 

r ' = 

f/(u o / L) 

N f.= N.*- N- + ./(m r L 2 ) 

r' = 

• • •• 2 ? 

*'= */ (UqL ) 

T '= T/(m r u 2 /L)(Tension) 

*9= 

Xg/L 

1 '= 1/L 

II 

V L 

x p= V L 

X A= 

x cr 

V (m r u o /L) 

X u/ m r 

yp= V L 
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